Thermal vortex dynamics in a two-dimensional condensate 
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We carry out an analytical and numerical study of the motion of an isolated vortex in thermal 
' equilibrium, the vortex being defined as the point singularity of a complex scalar field ■(/;(r, t) obeying 

I a nonlinear stochastic Schrodinger equation. Because hydrodynamic fluctuations are included in this 

description, the dynamical picture of the vortex emerges as that of both a massive particle in contact 
I with a heat bath, and as a passive scalar advected to a background random flow. We show that the 

vortex does not execute a simple random walk and that the probability distribution of vortex flights 
has non-Gaussian (exponential) tails. 
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I. INTRODUCTION 



CIh Vortices and other topological field configurations play a fundamental role in determining the properties of many 
: fascinating materials and control the physical mechanisms underlying several applications. Examples include su- 
perfluids 0, superconductors periodic solids liquid crystals |^, two-dimensional magnets (t), propagating 

coherent light beams |^ , and perhaps even the hot plasma that filled the very early Universe . Common aspects of 
the phenomenology of these diverse applications stem from the mathematical similarity of the underlying field theories 
^ ' and their classical static solutions. 
'■^ . Models of vortices as Brownian point particles, characterized by a mass, mutual interactions, and damping, are 
^ ■ commonly used to study the behavior of superconductors and superfluids . fn this picture the equation of motion 
O _ for a single overdamped vortex, not subject to any external forces, is extremely simple: 



7r(t) - f(t); 



{Mt)Mt')) = ^^s,,j{t-t'i (1) 



where 7 is the damping coefficient, M is the vortex mass, {(t) is a Gaussian thermal noise, fi^u G {x,y}, and (• • •) 
denotes canonical ensemble averaging. It follows that the vortex velocity distribution is also Gaussian. 
f-*) : This effective picture is extremely appealing — primarily because of the drastic reduction of the number of degrees 
of freedom. However, to our knowledge, it has not been shown to be derivable from the dynamics of an underlying 
^\ field theory. In this paper we study this question by solving both analytically and numerically for the motion of the 
vortex as an effective degree of freedom arising in a stochastic nonlinear Schrodinger equation. It is well-known that 
^ the conventional description of Brownian motion in terms of equations such as (0), that of a heavy particle interacting 
with light scatterers, ignores the presence of hydrodynamic fluctuations. A more complete description would allow us 
' not only to test the validity of the Brownian motion model for the vortex but also to compute corrections to it and, 
y possibly, to find new physical effects. 

Q Our analytic approach rests on the use of a singular perturbation expansion around a rigid vortex utilizing a low- 

^ temperature or weak-noise expansion. We are able to derive a Fokker-Planck equation for the single vortex distribution 
ILJ function which corresponds to that of a passive scalar advected to a background flow [not just a simple diffusion 
[ equation as would be predicted by Eqs. (|l])]. In our numerical work we are able to extract the diffusion constant 
for the vortex motion which turns out to be in good agreement with the theoretical prediction. Consistent with our 
?H ^ theoretical analysis, the vortex effective mass diverges logarithmically with the system size. We also investigate the 
statistics of vortex flights and demonstrate that the probability distribution function (PDF) has an exponential tail 
implying nontrivial correlations in the background thermal flow field. 

In Section I we present the stochastic nonlinear Schrodinger equation and discuss its physical relevance. In Section 
II we use the low-temperature perturbative expansion about the vortex collective coordinates to derive its equation 
of motion. Section III describes our numerical techniques and the associated results. We end with Section IV with a 
further discussion of our conclusions. 
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II. THE STOCHASTIC NONLINEAR SCHRODINGER EQUATION 



We study the motion of an isolated vortex by considering it as a singularity of a classical stochastic nonlinear 
Schrodinger field ■0(r,t) in two spatial dimensions. The equation of motion of this model theory is (a < 0): 

dib 

(7?(r, t)rj* (r', t')) = 2kBT^6{t - t')S{r - r'), 
(r?(r,tMr',O)-0. (2) 

The damping coefficient 7 and the stochastic force ?7(r,t) model the coupling of the condensate to a heat bath, 
satisfying the fluctuation-dissipation relation. Coefficients a, b and m are assumed constant and independent of 
temperature; thermal effects are fully described by the stochastic force. 

The above equations (|^) describe the phenomenology of physically important systems such as atomic Bose conden- 
sates and superfluid helium II [^,|3|. The same equation (at T = 0), amended as appropriate by the vector 
potential, describes the order-parameter dynamics of type-II superconductors . Although the above is known 
to be strictly true only for a class of type-II superconductors with magnetic impurities fl^ , it is often assumed that 
Eqs. (|^) apply to all type-II superconductors, unless indicated otherwise. We will henceforth study these equations 
as an important paradigmatic field theory representing a broad range of related physical phenomena. 

It is convenient to make all quantities dimcnsionless by standard substitutions. We set tp = 'ij)^—a/h, 7 = 7?i, 
t = i^/{—a), r = rh/\/—2am, 77 = fi{—a)yj—a/b, and ksT — —ah^ /{2j3mb). Note that r is normalized to a unit of 
length h/ \/—2am, which is the Ginzburg-Landau coherence length, usually denoted ^. With these substitutions, and 
after dropping the tildes, Eqs. (||) become 

P 

(77(r,i)ry(r',t'))=0. (3) 
The position of an isolated vortex, R(i) , is defined via a contour integral 

^[Varg-!/;(r,i)] • dl, (4) 

which equals 27r when the integration path encloses the point R(i), and vanishes otherwise. 

The Brownian motion of vortices as point particles has been argued to hold primarily in the overdamped limit 
(7 — > 00), which is directly relevant to the superconducting case; this is also the limit studied in this paper. We 
envision a situation in which a single vortex exists in the ground state of the field V'(r, t). In a thin superconducting 
film one can in principle achieve the same effect by placing the sample in a magnetic field of one flux quantum 
per sample area. In the latter case there are many other factors of practical importance, such as geometrical and 
point pinning, which however will not be dealt with here. The vortex is assumed to be located near the center 
of the spatial extent of the condensate, so that boundary effects can be neglected. The radius of the sample Ac 
provides a natural infrared cutoff. At sufficiently low temperature, thermally induced vortex-antivortex pairs will be 
exponentially suppressed and we may assume that there is only one vortex in a spatially bounded superfluid film at 
all times. This assumption has been verified numerically for a range of temperatures and sample sizes studied. 



III. PERTURBATIVE ANALYSIS OF VORTEX TRANSPORT 

In this Section we utilize the singular perturbation expansion of Kaup ||l^] to extract the vortex position from the 
stochastic equation of motion (^. We first expand the field as 

7/;(r,i) =^o(X) + e0(r,i); (5) 
^(r,i) =01 (X,i)+e(/.2(X, <) + •••, (6) 

where X = (AT, Y~) is the comoving coordinate, measured relative to a moving reference point; X(t) = r — R°(i). The 
absence of explicit time dependence in '!/'o(X) in the comoving reference frame indicates that it is the static vortex field 
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in the absence of fluctuations. Tomboulis has shown that introduction of a collective coordinate such as R°(i), 
while conserving the total number of dynamical degrees of freedom, is a canonical transformation. The fluctuating 
field (/)(r, t) can be thought of as a superposition of harmonic modes in the background of the rigid vortex ■00 (X). The 
phonons are gapless [[l9| , therefore they will be an important consideration at any finite temperature. The variable 
e ^ 1 plays the role of a small parameter as well as that of a bookkeeping device to control the perturbation series 
po| . Another useful way to think about e emerges when one absorbs e into the definition of ry(r, t) in Eq. (H): e then 
reappears as the square root of temperature in the fluctuation-dissipation relation in Eqs. (^). Therefore, the small-e 
expansion is actually the small temperature or weak noise expansion. The equation of motion (H) is required to hold 
at every order in e, which results in a hierarchy of equations for ipoCK), (f>i{X.,t), etc. 

Because of the added fluctuations in Eq. (||), the actual vortex position R(i) as the singularity of the full dynamical 
field V'lr, t) does not in general coincide with the coordinate R''(t), which by definition is the singularity of the rigidly 
moving static vortex field V'o(X). The two are related via 



Kit) = p[R°it), e<t)] = R°(t) + e / d^r 



<50(r) 



</)(r,t) + c.c. + 



(7) 



6=0 



where p[Il^{t),e<j>] is a vector fmictional of the fluctuation field e(f>{r,t) and of the reference point R''(i). We define 

R°(t) = + ev?(0+eMW + ---, (8) 



R(t) = + evi(t) + e2v2(t) 



(9) 



which makes explicit the fact that there is no vortex motion in the absence of thermal fluctuations. Differentiating 
Eq. with respect to t, we find at order e. 



vi(0 =v°(t)+V[R(0,t]. 



(10) 



where V[R(i),i] = J cPr "^^^"^(rj *^^ 4>{^i't) + cc. is a time-dependent velocity field. In the last expression we 

replaced R°(i) with R(i), as the two are equal to 0{e). The statistical properties of V[R(t), t] are given by those of 
the excitation field 0(r, i), i.e., the phonons, and will therefore depend on temperature. 

Substituting Eqs. (g), (^, (||) and into Eqs. (||) and collecting powers of e we obtain, at order e°, 

= -V^V^o-[l-|V'o|']^o. (11) 

The properties of the static vortex solution are well known; ipoO^) = /(|X|) exp[i aig{X + iY)]^ where f{x) ^ x for 
cc <C 1 and f{x) ~ 1 for x ^ 1. At order e we find 



d_ 

dt 



C 



Vx 



V'o 
V-o* 



(12) 



where £ is a Hermitean matrix 



L = 



*\2 



[l-2|Vo| 



■0^ 

[l-2|V^oP] 



(13) 



The noise 7?(X, t) remains white in space and time also in the moving reference frame, just as in the original Eqs. (|^). 
Differentiating Eq. (O) with respect to X, Y we find 



C 



dxro 



dy^O 



= 0. 



(14) 



Hence the null space of the linear operator C is spanned by the eigenvectors (dxipo,dx'4'o)'^ and (9y0Oi f^yV'o)^- 
These generate uniform translations of "00 in the xy plane and correspond to the Goldstone modes. 

It turns out |17| that all secular terms in the e-expansion vanish when we demand that the arbitrary perturbation 
(0, 0*)^ be orthogonal to the null space of the operator C. With this condition in place [which also uniquely specifies 
the reference point R°(i)] we now take the scalar product of Eq. ( p^ with the eigenvector {dx4'o,dx4>o)^ j 
integrate over the sample area to obtain 



= Re / d^X{dxro)Wi{t)-Vx^Jo-v{^,t)] 



(15) 
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The same equation is obtained if we use the other eigenvector in the scalar product, only now dx'4'o '^^ replaced by 
dY%. Solving Eq. (|l|) for v?(i) leads to 

v?(t)=F(0, (16) 
where the stochastic force F(i) is Gaussian with zero mean and autocorrelator 



{F,{t)F,{t')) = j^^^Kt ~ t')5,.. (17) 



where fi,v G {X, Y}, and 



Ac 



M(A,)= / d^X|Vx^o(X)|^ ~27rlnA, (18) 



can be interpreted as the (cutoff-dependent) inertial mass of the vortex. The fully dimensional form of this equation is 
(Ac) ~ (27r7V6)ln(Ac/C)- This quantity, in a different derivation, has received the same interpretation by Simanek 
The so-called vortex core mass |22| is implicitly present in the first part of Eq. (18) and makes a constant 
contribution to M(Ac), whereas the asymptotic expression in the latter part of Eq. (18) contains the Ac-dependence 
which is dominant in the limit of large Ac. In charged superfluids there is yet another contribution to the vortex 
inertial mass, the so-called electromagnetic mass ||2^ , which originates from the energy of the electric field generated 
by a moving vortex. This contribution is absent in Eq. (18), which strictly applies only to neutral superfluids. 

Substituting for v?(t) in Eq. (|l|) from Eq. (0) and using definitions (||) and (^, we arrive at an equation of motion 
for the vortex, valid to 0{e), 

jR{t)^eF{t) + eV[Rit),t]. (19) 

As the last step we set e = 1 in this equation. 

If we assume that V[R(t), i] is slowly varying relative to the noise F{t), then we may treat V[R(i), t] as an external 
field Q and arrive immediately at the Fokker-Planck equation |Q for the density c{r,t) = {6[r — R(i)]): 

+ V • [V(r, t)c(r, t)] = aV2c(r, t), (20) 

where the diffusion coefficient 

1 



/3M(Ac) ■ 



(21) 



This equation is of the same form as that for a passive scalar advected to background fluid flow V(r, t) P4| , p5| . It is 
well established — particularly in the best-studied case of incompressible fluid flow V • V(r, t) = — that in contrast to 
pure diffusion, an advected passive scalar may display non-Gaussian (for example, exponential) tails in the probability 
den sity function (PDF) c(r, t). These tails arise due to nontrivial spacetime correlations in the advecting velocity field 
1^,^. Interestingly, neither diffusion nor advection when acting alone can produce non-Gaussian tails in the PDF's. 
In our case, the non-Gaussian tails would appear in the PDF of vortex displacements during a fixed time interval At 
(vortex ffights). In order to test this possibility without detailed knowledge of the velocity field V(r, t), we study the 
thermal motion of an isolated vortex numerically as described below. 



IV. NUMERICAL ANALYSIS AND RESULTS 

In this section wc describe our numerical methods used to solve Eqs. (||) under the constraint that only a single 
vortex be present in the simulation volume. An unpaired vortex is forced into the sample by coupling the field 'ip{r,t) 
to a static external gauge potential A(r), replacing V in Eq. (6) by V-KiA(r). If we choose A(r) = — (27r/_L^)ya;, there 
will be on average one vortex per area L x L in an infinite 2D sample. The introduction of A(r) is a mathematical 
artifice and does not represent a dynamical field. In the language of superconductivity, A(r) represents an external 
uniform magnetic field that completely penetrates the sample without screening, and vanishes in the thermodynamic 
limit L oo. This is a realistic approximation to a superconducting thin film geometry, as in this case the external 
magnetic flux must fully penetrate the superconductor. Screening is only manifested as an inhomogeneity of the flux 
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density threading the plane of the superconductor. The characteristic length scale for these variations is l|2^ A = A^/d, 
where A is the bulk magnetic penetration length of the superconducting material, and d is the thickness of the film. 
We study the system in the limit where A is much larger than any other length of the problem. In this limit the 
flux density threading the superconductor is uniform and screening can be neglected. One might also worry that the 
derivation leading up to Eq. (18), which was done explicitly for the case of the neutral superfluid, no longer applies 
when A(r) is introduced. While it is a fact that the massless Goldstone boson in a superluid is in general replaced by 
a massive vector boson in a superconductor, the mass of this boson vanishes in the thin film limit considered, and the 
gapless phonon spectrum is restored, as pointed out by Simanek. [2l| ] By the same token the electromagnetic vortex 
mass vanishes in our thin film limit, as the moving vortex does not generate any electric field. 

We choose our system to be a square, S = (0,1/) x (0,1/). We wish to employ a variant of periodic boundary 
conditions so that the entire xy plane consists of physically equivalent copies of the system at all times; by this we 
mean that the local superfluid density |i/)(r, i)p, and the supercurrent density 

J(r,i) = ^(^VV*-V'*VV') + |V'l'A (22) 

should both be periodic functions with periodicity L in both spatial directions. To this end we impose quasi-periodic 
boundary conditions on ip{r,t) (r S S, Uy and riy arbitrary integers): 

'0(r + rixLx + UyLy, t) = ■(/'(r, t) explin^Ox + iny{2TTx/ L + Oy)]. (23) 

The angles 9^ and 6y are as yet arbitrary and may be functions of time. Upon discretization of Eqs. (H) we observe 
that any given time dependence of 9x and 9y will in general affect the motion of ■0(r,t). Yet 9x and 9y are not 
dynamical degrees of freedom in the sense that their motion is not determined by Eqs. (H). Therefore, we are left with 
a problem of choosing the "right" time dependence for 9^ and 9y. Fortunately the physical interpretation of these 
variables is straightforward. In order to find it, we rewrite the stochastic field equation in (|^) in the form 

(V7 - 1)^^ = ^"^^J^ ^ ^ + «/(r, t), (24) 

where 

o ...o 1 



-|(V + ^A)^r-|V^|^ + -|^|4 



(25) 



is the Hamiltonian of the system (the condensate) disconnected from the thermal reservoir. One can show that 

= x - [ dy J.(r,t)U^o. (26) 
1 Jo 



dH „ r SH 

The first part of this equation is an algebraic equality; in the second part we used the definition of the supercurrent, 
whose explicit form appears in Eq. (|2j). Therefore, the first derivative of the Hamiltonian with respect to 9x is equal 
to the net supercurrent passing through the edge of the sample perpendicular to the x axis. A completely analogous 
statement holds for the other spatial direction as well. 

If one is only interested in the thermal motion of the vortex in the absence of external forces — such as, for example, 
the Magnus force created by macroscopic superflow — the desired constraints on angles 9x and 9y are 

These constraints define a zero flow ensemble p7| , relevant for experiments with no externally imposed driving flow. 
The angles 9x and 9y, though not dynamical variables, have their values fixed by two dynamical constraints, (|27|), at 
all times. 

As a result of the quasi-periodic boundary conditions, the isolated vortex is part of an infinite (always perfectly 
square) vortex lattice, with lattice constant L. The physical cutoff A^ needed to evaluate the vortex mass is estimated 
as Ac ~ L/ y/iT (this is when a square of side L and a circle of radius Ac have equal areas). Our numerical approach 
consists of simultaneously solving the lattice version of Eqs. (^ on the spatial domain 5, and the two dynamical con- 
straints (p7|). The numerical grid spacing 5 must resolve both the vortex core size f , which is the unit of length, and the 
phonon correlation length ^p, which is the characteristic length scale of the four-point correlator (|'0(r, <)p|^/'(r', t)p), 
and depends on temperature. We have found it adequate to set 5 = 0.15 throughout the temperature range studied. 
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Numerical stability of the second-order stochastic Runge-Kutta method |28| used to integrate Eqs. (|^) then allows a 
time step of r = 0.005. We study sample sizes L = 6, 12, and 24, and temperatures 1//3 = 0.0125, 0.025, 0.05, and 
0.1. 

The vortex is clearly identifiable and no thermal vortex- antivortex pairs are generated during the numerical runs. 
These ranged in duration from 6 x 10"* to 6 x 10^, in dimensionless units (up to 1.2 x 10^ time steps). The characteristic 
equilibration time was ^ 500. The vortex position is taken to be the coordinate of the lower left-hand corner of the 
plaquette where it is detected with the discretized version of the integral (|^). 

Fig. 1 shows ([AR(i)]^) for three different L and the same (3. It is apparent from the data offsets that vortex 
motion is not purely diffusive. On the other hand, the calculated diffusion coefficient a correctly approximates the 
slopes of the numerical data at long times. Therefore we conclude that the vortex mass, measured through cr, displays 
a logarithmic divergence with the system size L. The issue of the value of the vortex inertial mass has been the 
subject of much controversy in the literature [ p9|j30| . Our results show that even in the presence of random flow the 
logarithmic divergence of the classical vortex mass persists. This stands in agreement with an earlier result in 
the presence of coherent flow. It would be interesting to investigate whether this divergence persists in the presence 
of multiple vortices or in the limit of low damping. 

An independent measure of diffusion is obtained through the concept of the mean first passage time t(£) over the 
boundaries of the square box of size 2t x 21. For pure diffusion one has = 0.295 1"^ /a Fig. 2 reiterates that the 
vortex motion is diffusive at long times. Deviations from pure diffusion are manifest at short times that correspond 
to vortex flights of 0(1) and less. The evidence from Fig. 2 also shows that a depends linearly on temperature, as 
expected from Eq. (21). This is a useful check on the quality of our numerics — the vortex is not subject to spurious 
pinning due to the discreteness of the underlying lattice. The flattening-out of the data for large £ is due to finite 
sampling. 

Finally in Fig. 3 we present the PDF of vortex flights. This distribution has exponential tails (straight lines on 
a logarithmic scale as in Fig. 3), in addition to the expected Gaussian profile typical of Brownian motion. Our 
numerical results show that these tails persist for all system sizes and (3 studied, and the slope decreases as the vortex 
mass A/(Ac) is increased (which corresponds to a volume increase) . The presence of exponential tails has been linked 
to nontrivial time correlations in the advecting velocity field p4|. The presence of finite correlation times can be 
either a feature of a stochastic velocity field or can arise from deterministic, chaotic motion (e.g., intermittency). In 
our case, we expect the exponential tails to arise from the time correlations of overdamped phonon modes which are 
allowed to exist in the fluctuation spectrum of our theory. This aspect is presently under investigation. 



V. DISCUSSION AND CONCLUSIONS 



The particular models of vortex motion of the type described in the Introduction typically ignore the presence of 
hydrodynamic modes. In contrast, real physical systems, in particular those described by stochastic field theories, 
necessarily include these effects. Consequently, as demonstrated here both analytically and numerically, the motion of 
a single vortex must be viewed as diffusive and randomly advected at the same time. The Fokker-Planck equation ( [20| ) 
makes explicit the formal analogy between single vortex fluctuations and the dynamics of a passive scalar advected 
to random fluid flow. The quantity c(r, t) is the ensemble-averaged probability density of finding a vortex at r at an 
instant t. The ensemble averaging is done over the fast degrees of freedom of the heat bath. The slow, hydrodynamic, 
degrees of freedom of the superfluid — the phonons — remain implicitly present in the velocity field V(r,<:). This 
equation presents a picture of the vortex as a diffusing massive particle subject to random advection. The advection 
is apparently due to phonons, which propagate through the vortex core randomly. As a result, the probability of 
large-scale vortex flights is increased. 

In turbulent fluid flow experiments with a dye as the passive scalar it is possible to measure the dye concentration, 
analogous to our c(r, t), directly, because individual dye particles do not interact among themselves (apart from 
contact interaction). Vortices, however, interact logarithmically, so one cannot conflrm the exponential tails in c(r, t) 
directly by studying motion of a large aggregate of vortices in a sample. Instead, one must study the motion of a 
single vortex, and measure the statistics of vortex flights. 

Examples where thermal vortex motion plays a key role include theories of thermal depinning, flux creep, and dilute 
vortex lattice melting in high- Tc superconductors. Our current understanding of these phenomena relies exclusively 
on the Brownian particle description. It would be interesting to re-examine these problems in the light of the 
advected scalar aspect of vortex motion. An exciting new possibility is a direct experimental observation of vortex 
motion in two-dimensional high- Tc superconductors by fast scanning tunneling microscopy Js^ . With this method 
the overwhelmingly more frequent occurrence of large vortex flights in comparison with that of a simple Brownian 
particle could be put to direct test. 
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100 200 300 400 

t 

FIG. 1. The mean-squared vortex displacement vs. time, for three different sample sizes (top to bottom): L — 6, 12, and 
24. The temperature 1//3 — 0.1 is the same for all plots. The straight lines are the Einstein formula for pure diffusion, 
([R(t) - R(0)]^> = 4at, with a calculated from Eq. (|ip. 




FIG. 2. The mean first passage time vs. the box size, for four temperatures (top to bottom): 1//3 = 0.0125, 0.025, 0.05, 
and 0.1. System size L = 6 is the same for all plots. Straight lines are the theoretical predictions for pure diffusion, with a 
calculated from Eq. (pl[). 
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FIG. 3. The PDF of vortex flights, for L = 6 (triangles) and L — 12 (squares, shifted by four orders for clarity). = 0.1, 
and time interval between measurements At = 0.1. Full and empty symbols correspond to V{Ax) and ViAy). The dotted line 
is the Gaussian distribution with the same rms deviation as the numerical data for L — 6. The straight lines are drawn with 
slopes corresponding to numerical data. 
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